Dynamical polarizability of graphene irradiated by circularly polarized ac electric 
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We examine the low-energy physics of graphene in the presence of a circularly polarized elec- 
tric field in the terahertz regime. Specifically, we derive a general expression for the dynamical 
polarizability of graphene irradiated by an ac electric field. Several approximations are developed 
that allow one to develop a semianalytical theory for the weak field regime. The ac field changes 
qualitatively the single and many electron excitations of graphene: undoped samples may exhibit 
collective excitations (in contrast to the equilibrium situation), and the properties of the excitations 
in doped graphene are strongly influenced by the ac field. We also show that the intensity of the 
external field is the critical control parameter for the stability of these excitations. 

PACS numbers: 



I. INTRODUCTION 

Graphene is a genuinely two dimensional material 
whose peculiar properties have received a lot of atten- 
tion since its first isolation in 2004. 1,2 Structurally, it 
is a single atom thick layer of graphite, i.e. a two di- 
mensional crystal, that remains stable both when it is 
deposited over a substrate or when it is suspended. Its 
electronic properties have attracted huge interest: the 
low energy excitations are chiral massless Dirac elec- 
trons in two dimensions, thereby providing a new plat- 
form for testing the basic tenets of solid state physics. 
This fact, which ultimately arises from the honeycomb 
structure of the graphene crystal lattice, is responsible 
for a strikingly different electronic behavior as compared 
to the conventional two dimensional electron gases (e.g., 
in semiconductor heterostructures) studied extensively in 
the laboratory.^! 

The effect of external fields in the low energy proper- 
ties of the electric carriers in graphene has been a topic 
of extensive research since early days, as the discovery 
or the anomalous Quantum Hall effect witnesses Un- 
derstanding the behavior of graphene in the presence of 
electrical and magnetic fields is of major relevance both 
from a fundamental and an applied point of view. The 
former, since new exotic behavior may arise in the pres- 
ence of external fields, and the latter, because external 
fields can be used to manipulate its properties, for in- 
stance by opening gaps in the electronic spectrum, which 
is essential for applications in the semiconductor indus- 
try. 

The effect of radiation on both monolayer and bilayer 
graphene has been analyzed only recently, and has led 
to the prediction of a variety of phenomena, such as 
the photovoltaic Hall effect,— metal-insulator transition 
of graphene^ valley-polarized currents in both monolayer 
and bilayer graphene^-'- and photoinduced quantum Hall 



effect in the absence of magnetic fields^ Other theo- 
retical works include the analysis of ac transport prop- 
erties through graphene ribbonSj 11 graphene-based in- 
junctions^ graphene-based Fabry-Perot devices^ and 
the recent proposal of quantum pumping in graphene by 
an external ac field. 14 Experimentally it has been found 
that a circularly polarized ac field induces a dynamic Hall 
effect in graphene.— Several studies have been devoted 
to the theoretical analysis of the quasienergy spectrum 
of graphene and graphene dots under ac fields ^r— and 
the optical properties of graphene have been studied by 
calculating the optical conductivity^ One of the earli- 
est and yet most important findings in all these studies 
is that a circularly polarized field induces a band-gap 
at the Dirac point, along with dynamical gaps at other 
momenta, all of which are tunable by the field intensity. 
This is, however, not the case for a linearly polarized 
field: there the anisotropic quasienergy spectrum shows 
dynamical gaps at non-zero momentum only in certain 
directions, and especially no gap is induced at the Dirac 
pointM 

In this paper we study theoretically the effect of a cir- 
cularly polarized ac electric field in the terahertz regime 
on the electron excitation spectrum, and on the electron- 
electron interaction. The interactions are found to be 
affected qualitatively by the external field, altering the 
nature of the single particle excitations as well as the 
many-particle excitations, both in doped and undoped 
graphene. Special attention is paid to the existence of 
a plasmon in undoped graphene, which is not present in 
its field-free counterpart. In order to perform this in- 
vestigation, the natural object to study is the dynamical 
polarizability, which has already been studied extensively 
in graphene without an ac field<2i _ — 

The structure of the paper is the following: In Sect. HT1 
we briefly introduce the Hamiltonian of graphene in the 
presence of a circularly polarized electric field, emphasiz- 
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ing the role of Floquet theory in Sect. Ill A[ and present 
several approximations to the single electron Hamilto- 
nian valid for weak fields in Sect. Ill Bl Section I1III is 
dedicated to the analysis of the dynamical polarizability: 
we derive a general expression for the polarizability of 
graphene in an ac electric field in Sect. lIII Al and compare 
it with the corresponding expression for the two dimen- 
sional electron gas^ Finally in section UlIBI the general 
formula is combined with the analytical approximations 
in order to work it out both for the non-interacting sys- 
tem and for the interacting system in the Random Phase 
Approximation (RPA). 

II. SINGLE ELECTRON PROPERTIES OF 
GRAPHENE UNDER A CIRCULARLY 
POLARIZED AC ELECTRIC FIELD 

A. Model and technique 

In the low-energy regime, the Hamiltonian for single 
electron excitations in graphene is the infamous Dirac 
Hamiltonian. In order to introduce a time-dependent 
electric field we choose a gauge in which the latter is 
represented via a gauge potential A(t), whose time de- 
pendence is that of a single monochromatic and circularly 
polarized wave of frequency O: 

A(t) - -^=[xsm(Qt) -ycos(Ot)], (1) 

By using a minimal coupling scheme, the Hamiltonian 
for graphene irradiated by this electric field reads: 

_ / k x - ik y + iAe- int \ 

H[t) ~ [k x +ik y -iAe mt ) ' {1 > 

with A = eE /{s/2Q) and uf = H = 1. Here the 
Hamiltonian is expressed in terms of Bloch states of 
momentum k, which is defined with respect to one of 
the valleys. Note that the electric field does not couple 
the spin and valley degrees of freedom in graphene, 
which remain as an extra degeneracy N V N S — 4. The 
eigenstates of graphene in the absence of the external 
field are two-dimensional spinors representing the two 
components of the unit cell of the honeycomb lattice 
in graphene, that once diagonalized give rise to two 
bands (or Dirac cones). In the Dirac Hamiltonian, the 
pseudospin has a scalar coupling with the momentum, 
and its eigenstates are those whose pseudospin is either 
parallel or antiparallel to its momentum. In fact, the 
mathematical structure of the Hamiltonian coincides 
with that of an electronic spin coupled through Rashba 
interaction to a magnetic field. In this analogy, the 
momentum in graphene plays the role of the magnetic 
field, and the pseudospin operator is the analogous to 
the ordinary spin, both having the same representation 
in terms of Pauli matrices. This allows one to write 
the Hamiltonian H = a • k. In the presence of an 



external electric field an extra term of the same nature 
arises in the Hamiltonian, now coupling the pseudospin 
and the electric field, and inducing transitions between 
the eigenstates for the isolated system. In a sense 
the momentum and the electric field are competing 
dynamically for the direction of the pseudospin, but no 
compromise can be reached due to the time-dependence 
of the field, which no longer allows for an analysis of the 
problem in terms of stationary eigenstates. 

To proceed, we apply the Floquet theorem, which is 
the most suitable way to address time periodic Hamilto- 
nians (detailed accounts can be found in Refs. [29l - [3lj | ) . 
Floquet theory states that for a Hamiltonian that is pe- 
riodic in time - H(t + 27r/f2) = H(t) - a complete set of 
solutions of the time-dependent Schrodinger equation 

H{t)m))=ij t \m) (3) 

can be written as 

\cf> a (t)) = \<t>a(t + T)), (4) 

where a contains the quantum numbers of the problem 
and the so-called Floquet index, that we will label as I. 
The role of this index is to classify the different sidebands, 
since e Q , the quasienergies, are defined mod Ml, being 
related by the simple transformation: 

e Q (z) = e Q (o) + '0 (5) 

In analogy to the Bloch theorem, the quasienergies can 
be mapped into a first time Brillouin zone, which is 
[— 0/2, 0/2], and therefore corresponds to I = 0. 

The Floquet states \<fi a (t)} have the same periodicity 
as the driving field (see Eq.(j4])) and can therefore be ex- 
panded into a Fourier series: 

oo 

\Mt))= £ e' ; " ot |C>- (6) 

n— — oc 

The Floquet states are also defined in different branches 
of solutions, being related between them by the transfor- 
mation: 

l€(o> = l*S)> (?) 

Substituting Eq.(|3]) into Eq.([3]) and using Eq.© yields a 
static eigenvalue equation of the form 

Y l (H nm -nnS mn )\4,^)=e a \<f>l). (8) 

m 

Defining now the Floquet Hamiltonian as 
Hp m = H nm — mtt5 mn , we see that a significant 
simplification has been achieved: the time-dependent 
problem has been transformed to a static problem, and, 



3 



consequently, one can apply the intuition about equilib- 
rium problems to make statements about a dynamical 
problem. The resulting equilibrium-like observables 
derived within this framework have to be understood as 
time averages over a period of the external field. 

Let us now apply the Floquet formalism to the Hamil- 
tonian of graphene @. In this case, the solutions are 
characterized by indices a = (k, a, I), being a = ± the 
pseudospin index: 

nn^ a + (k x - ik y )cj>2 b + iA0 n a +hb = e Q C a 
nf>e" + (*» + ihy)4>T ~ iMl~ 1,a = £aC' b (9) 



Notice that a and b are the indices for the sublattices of 
the honeycomb lattice. These equations can be written 
in matrix form, where the infinite Floquet Hamiltonian 
reads 
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(10) 



The structure of this Hamiltonian deserves a few com- 
ments. The ac field A connects (2 x 2) graphene Hamil- 
tonians with energies nQ, and (n+l)fl and so on. Each of 
these building blocks contributes with its own dispersion 
relation (that of a Dirac cone) to the energy spectrum, 
and the field introduces transitions between these cones. 
These transitions are expressed as anticrossings in the 
spectrum, and become exact crossings for A — > 0. It can 
easily be seen that the anticrossings occur at |k| w nil/2, 
with n = 0,1,2,.... At |k| w 0/2 e.g., the (+,0) and 
the (—,1) sideband anticross, which would be a so-called 
one-photon resonance. 
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B. Analytical approximations to the single particle 
Hamiltonian 



The Floquet Hamiltonian, Eq. (ITU|) can be diagonalized 
numerically in order to analyze the energy spectrum and 
its features. However, in order to simplify calculations 
and to illuminate the main physics, here we resort to an- 
alytical approximations which capture the main features 
whenever the electric field intensity is sufficiently weak. 
We will show in Sect. Ill CI the full numerical results for 
the quasienergy spectrum in order to compare it with 
the analytical approximations. Before introducing such 
approximations, it is convenient to project the Hamilto- 
nian (|TU)) into another basis. As can be seen in Fig. [TJ 



FIG. 1: Sketch of the Hamiltonian for the circularly polarized 
field. Note that if k = 0, the Hamiltonian breaks up into 
disconnected two-level systems, in which site a„ is coupled to 
site fe„+i. 



for k = 0, the Floquet chain breaks up into a series of 
disconnected two-level systems. We therefore diagonal- 
ize the Hamiltonian for the two-level-system, and then 
write the full Hamiltonian in the resulting basis. The 
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Hamiltonian for k = reads 



^o)(CL a ol + l^o)(CL b l 



An excerpt of the series of (2 x 2) Hamiltonians is 



(11) 



rrk=0 

tip 





iA 








-iA 


(n + i)n 














(n-l)Cl 


iA 








-iA 





(12) 



Out of the four eigenenergies of this matrix we are inter- 
ested in 



with 



e± = ICl ± - A, 



n = ^aa 2 + n 2 . 



(13) 



These two energies fulfill liniyi^o = ^, thus, we as- 
sociate the first Brillouin zone, I = 0, from the Floquet 
solutions, to the solutions corresponding to graphene in 
the absence of an external field. The corresponding eigen- 
vectors are 



(a + 2O)|^L ; 



\<t>D = ^ ((A + 2fi)|^ > + 2iA\^ b ) 



(14) 
(15) 



where N = \J AA 2 + (A + 2fl) 2 . From here on and for 
the rest of the paper, we neglect the index k in the ener- 
gies and vectors, unless we have to distinguish between 
k and k + q. 

Using these eigenvectors as a basis, we rewrite the full 
Floquet Hamiltonian (fTU|) 
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(16) 



r 



where fc = |k|, O = aTctanky/k x and we introduced the 
three functions Fq , F\ and F 2 that will form the basis for 
our approximations: 



Fx = 



(A + 2n) 2 



Fo 



AA 2 + (A - 


h2ft) 2 


2iA(A + 


2Q) 


AA 2 + (A - 


h2ft) 2 


4A 2 




AA 2 + (A - 


h2ft) 2 



(17) 



All three functions in Eq. (fTT|) are functions of A and f2. 
For small A/Q <C 1, one finds that F « 1, and Fi 2 ~ 
0, however, F± increases linearly whereas F 2 increases 



quadratically with A. Note that for a two-level system 
driven by a linearly polarized field, the nth order Bessel 
function J n (A/ft) plays the role of the function F ^ :2 
presented here, and for a complete analysis one has to 
consider Bessel functions up to infinite order, see e.g. 
Ref. [Uli. Here however, the complete information lies 
in -Fb,i.2- In the subsequent analysis, we will at first only 
consider the couplings given by F , and then include also 
the couplings given by F\ . We will neglect F 2 in general, 
which is valid for small A/Q,. 
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1. Fo-approximation 

A first approximation consists in neglecting both Fx 
and F 2 and considering only F 0l which connects energies 
with the same photon number n. This approximation is 
valid for the calculation of many observables as far as the 
dimensionless quantity A/Q -C 1 - i.e. the field intensity 
is small compared to the frequency - and we are inter- 
ested in excitations in the low energy sector, as we will 
see below, when we analyze the excitation spectrum and 
the generalized density of states. The resulting Hamil- 
tonian (fT6|) is then block diagonal with building blocks 



by varying the field strength of the applied ac field, 
see also Refs. [11320]. We point out that an analo- 
gous phenomenon occurs in the optics of semiconduc- 
tors in a strong THz-field: there the dynamical Franz- 
Keldysh effec t 33 ' 34 blue-shifts the conduction band edge 
(or, equivalently, the optical absorption edge) by the pon- 
deromotive energy, which also depends quadratically on 
the ac-field amplitude. This so called Fo-approximation 
neglects the coupling between Hamiltonians with a dif- 
ferent number of photons, and is therefore not useful 
once we are interested in the anticrossings of the Floquet 
quasienergies for non-zero momentum. 



HS , where the matrix HS reads 

t ' CO 



r 



4 W e 

F ke~* e e- , 



Its eigenvalues and eigenvectors are 



1 



e± Fo = ZO±-^4F 2 fc 2 + A 2 



(18) 



(19) 



2. F\ -approximation 

In order to analyze higher order processes, we go one 
step further and take into account the coupling elements 
F\, which capture the one-photon resonances, yielding 
a much more robust approximation for the Hamiltonian 
H-p (Unj) • At the resonances the relevant couplings are the 



ones between e„ —1 and e 



t , e+ and e n+1 etc., see Eq. (fl"6|) 



By applying the unitary matrix that diagonalizes Hp 



where 



V\xa\ 



Xb 



\Xb\ 



{Xa\4>t)+Xb\4>i)) 

{xt\4>t) - xl\<t>7)) > (20) 



Xa = 2F ke^ 
Xb = 



4F 2 k 2 + A 2 



(21) 
(22) 



The main virtue of this approximation is the fact that it 
captures the gap A produced at k = by the ac elec- 
tric field, giving an analytical expression for its magni- 
tude, A = y/4A 2 + O 2 - O, so this gap can be tuned 



Un 



v\xa\ 



\Xb\ 



Xa X b 
,Xb -Xa, 



(23) 



we can construct a new effective infinite Hamiltonian 
which includes the features of the one-photon resonance, 
and which is again block diagonal, now mixing the sectors 
that differ in one photon in the ^"approximation: 



7-7-cff.n 



fc n-l,F 

2 F F*k 2 e- 



-^FoFik-e 



(24) 



where S& = ^/4F^k 2 + A 2 . The intensity of the coupling 
is proportional to Fx, as expected. Diagonalizing this 
Hamiltonian yields the following Floquet quasienergies: 



in + i (o - y J(n-s k ) 2 + ^F3\F 1 \*k* > ) ak<k c 
(i + 1)0 - ± (o - J (n - s k f + §F2\Fx\*kA iffc>fc c 



in - \ o - J{n - s k y + mF 2 \Fx\ 2 k 4 ) ak<k 



{i - 1)0 + \ (o - Mn - s k y + i|F 2 |Fi| 2 fc 4 if k > k c , 



(26) 

where k c is the momentum at which the one-photon resonance takes place. For the Floquet eigenvectors, it is more 
convenient to write them in the basis that diagonalizes the Hamiltonian for k = 0, reading 

F 1 Tm^r^hmm & Xa ^ +%xmt) +cxt\<t>t +1 ) -&6i*r +l >] if k > ^ 
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|£- ) = J V(IXa| 2 + |Xb| 2 )(|«a| 2 + |«b| 2 ) 
VdXaP + IXbl^d^P + kbl 2 ) 



[eaX a \4ti) + €°Xb\4>T-i) -f6Xil^) + 6x:i0r)] if fc > *« 



(28) 



where we have introduced 




(n - s fc ) 



16 



F 2 |^| 2 fc 4 e*' 



6 = ^Fbi^e"^. 



(29) 
(30) 



The F\ -approximation captures the gap at k = as well 
as the first resonance. The latter gives rise to the opening 
of new gaps, whose expression can be obtained analyti- 
cally in this approximation, yielding 



16 



a Fi = J(s kc -ny + w F^F?ki 



(31) 



For a frequency ~ 150meV in the mid- infrared regime, 
and field intensity E « 4.8MV/m (so that A/Q = 0.1), 
the size of the two gaps would be A s» 3meV, A Fl » 
15meV. 



-H 



-H 



1 



-1 
1 



-1 
1 



-1 




r 





k r 



C. Single particle properties of the Hamiltonian 
derived from the analytical approximations 

We next consider the quasienergy spectrum for circu- 
larly polarized field, both the full numerical result and 
the analytical approximation derived in the previous sec- 
tions. With increasing field strength, zero-photon, one- 
photon, two-photon and higher order resonances appear. 
In Fig. [5] we compare the numerical (upper panel) and 
the analytical results (middle and lower panels) for the 
quasienergy spectrum as a function of the wavevector k x 
and for weak fields. In the middle panel we plot the 
-Fb-approximation, which reproduces very well the gap 
at k x = 0, but no other features induced by the ac field 
show up. The lower panel shows the results for the F\- 
approximation, which in addition to the Fo result cap- 
tures nicely the one-photon resonance at k c ps ±0.5. 

The analytical approximations can be tested by com- 
puting the density of states (DOS), which was already 
analyzed numerically using the full Floquct Hamiltonian 
by Oka et al.fi- Calvo et al., 16 and Zhou et al.M The 
generalized (time averaged) density of states can be cal- 
culated as: 

D(w) = 4^%-e k)CT ,o) ) (32) 

k,(r 

where the quasienergies are those defined in the first Bril- 
louin zone of the Floquet spectrum. In the -Fo - approxi- 
mation, ek,o-,z=o = e^Fo' see Eq. (fl9)) . and the density of 



FIG. 2: Quasienergy spectrum as a function of k x for k y = 0. 
The solid lines represent the / = band, the dashed lines the 
1 = 1, and the dashdotted lines the I = — 1 sideband. Upper 
panel: The full numerical result of the quasienergies. Middle 
panel: The quasienergies for the zero-photon approximation 
Fo. Lower panel: The quasienergies for the one-photon ap- 
proximation Fi. Parameters: k y — 0, A = 0.2, Q — 1. 

states can be calculated analytically yielding 

£ Fo M = ^Me(V|-^. (33) 

Notice the presence of the gap at zero energy in the den- 
sity of states. 

The simplicity of this zero-photon approximation al- 
lows for analytical computations of many physical quan- 
tities, something that no longer happens in general in 
the one-photon approximation, for which we have to re- 
sort to numerical calculations in most of the cases. For 
the generalized density of states, by using the analyti- 
cal quasienergies, Eqs. (|25p and (J2H]), the results of both 
the Fq- and Fi-approximation are plotted in Fig. [3] (up- 
per panel) . As a comparison, the lower panel of Fig. [3] 
shows the density of states calculated numerically by 
diagonalizing the full Floquet Hamiltonian (fTD|) . Once 
again, notice that the Fn-approximation works very well 
for energies of the order of the first gap, while in order 
to study the first resonance the F\ -approximation excels 
quite well. Higher resonances - visible in the numerical 
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result for the density of states at around u> w 1, are al- 
most negligible for the field strength we are considering 
here. 

Although the density of states has already been achieved 
numerically by various authors ^ 16 ' 20 the developed ap- 
proximations are useful in order to both obtain analytical 
results when this is possible or at least simplify the nu- 
merical complexity of the problem, as it happens when 
we use the one-photon approximations. In principle, 
these approximations are valid for arbitrary k, as long 
as A /CI -c 1. The results for the density of states show 
that in both analytical approximations and in the full 
numerical case the gaps remain stable independently of 
the range of integration in the momentum included in the 
density of states, i.e., the inclusion of higher momentum 
states does not close th e g aps in our case, contrary to 
what was found in Ref. [201 ] . 
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FIG. 3: Density of states versus energy for the analytical 
approximations Fq and Fi (upper panel) and considering the 
full Floquet Hamiltonian Eq. (fl0j) (lower panel). The Fo- 
approximation reproduces the gap given by A, see text. For 
the one-photon resonance in the Fi-approximation, the gap 
at ui — 0.5 is reproduced, where the field couples modes with 
n and n + 1 photons. For the same field strength as considered 
in the Fq and Fi approximations, the full numerical density 
of states is identical up to an additional resonance at oj ~ 
1, which is due to two-photon processes. In the inset, the 
region around the gap ui — is blown up for better visibility. 
Parameters: A = 0.1, fi = 1. 



III. SINGLE AND MANY-PARTICLE 
EXCITATIONS IN GRAPHENE IN A 
CIRCULARLY POLARIZED AC ELECTRIC 
FIELD 

A. Electron interactions and the formula for the 
dynamical polarizability 

So far we have analyzed the single particle properties 
of the Hamiltonian. However, a full description of elec- 
tron excitations in graphene requires to understand the 
role of electron-electron interactions in the system. The 
Hamiltonian of the interacting system in the presence of 



an ac field reads now, in second quantization, 
H(t) =v F y ^<t • (k - eA(t))¥ k 



k 



w q n q n, 



(34) 



where u q = 2ne 2 /eoq is the 2D unscreened Coulomb 
interaction. In the absence of the external ac field, 
this Hamiltonian has been extensively studied (for a re- 
view, see ref. |H). For doped samples of graphene, the 
Coulomb interaction becomes screened, yielding a system 
whose low-energy excitations around the Fermi surface 
are barely interacting, i.e., electrons in graphene behave 
as a Fermi liquid. Moreover, a collective excitation, a 
plasmon, exists j 22 i 23 This is no longer true when the level 
of doping is zero or very small, where the role of inter- 
actions is controversial due to the singular nature of the 
Dirac point, where screening is uneffectivej 24 ' 36 In order 
to understand the effect of interactions between electrons 
when an external ac field is applied, we will compute 
the dynamical polarizability, which tells us about the re- 
sponse of the system to probes that couple to the electric 
charge. This function will yield information about the 
full spectrum of electron interactions of the system, that 
contains both single particle and collective excitations. 
The dynamical polarizability of graphene in the presence 
of an ac electric field shows particular features that differ 
from its counterpart, the two dimensional electron g as, as 
well as from the one derived for graphene in the absence 
of the fieldi^i The detailed derivation of the polarizabil- 
ity function is given in the Appendix, whereas here we 
present the final result: 



a, a' k I 



f\x.,(T /k+q,( 



6 k.a,0 



/ n,b,* A n ib 1 2 



(35) 



The index n stands for the Fourier component of a solu- 
tion in the sideband I of the infinite Floquet Hamiltonian. 
The summation over n constitutes the scalar product of 
the solution \<ftk+q,cr' ,i) with |</>k,er.o), where we have used 
the fact that solutions belonging to the Ith Brilloun zone 
in the Floquet spectrum are those of the first Brilloun 
zone shifted by I units, see Eqs. §5§ and ([7]). The solu- 
tion I = is the one which fulfills the condition that at 
A — >• 0, n(q, ui) becomes the polarizability for an isolated 
graphene sheet. Mathematically, it is important to notice 
that the analytical properties of the dynamical polariz- 
ability do not change in the presence of the external ac 
field. It is a complex function, analytical in the upper 
half plane, whose real and imaginary parts are not inde- 
pendent, but related via the Kramers-Kronig relations, 
see e.g. ref. [37[ The latter can also be seen as a conse- 
quence of the causality in the response of the system to 
the external probe. 

The polarizability is written in terms of the single par- 
ticle excitations of the system, as in conventional linear 
response theory. In the presence of the ac electric field, 
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there is an infinite set of single particle excitations that 
differ between them in the relative number of photons. 
Once the external field is switched on, the system is no 
longer isolated, and the field can pump or extract energy 
into the system in the form of photons of frequency fi. 
Therefore, the polarizability can be seen as a linear com- 
bination of polarizabilities, each describing excitations in 
which the number of photons in the system changes by 
a certain integer number I. Or in other words, the re- 
sponse of the system to an external probe of energy u) 
and momentum q can arise from excitations in which no 
extra photons are introduced in the system, as it hap- 
pens in the absence of the field, or also in which a given 
number I of photons is introduced or extracted from it. 
Similar results as those shown here have been derived 
in the context of low-dimensional semiconductors^ and 
the 2DEG. 2 - The latter is usually the benchmark to com- 
pare the results derived for graphene, and therefore we 
include here the formula for the ease of comparison: 



n 



2DEG 



/k — /k+c 



k I 



W - £k+q,Z + £k,0 + IT] 



(36) 



Note in this expression that the effects due to the ac 
electric field appear in two different places: (i) the index 
I of the sideband reflecting a change in the number of 
photons, and (ii) as a modification of the single particle 
excitations in ek, CT .z- The situation is more complicated 
in the case of graphene, Eq. (|35|) . where also a momen- 
tum dependent overlap term between the excitations with 
different momentum must be included. This overlap is 
reflected in the polarizability via the scalar product be- 
tween quasieigenstates, and it is in turn a consequence of 
the existence of the pseudospin in graphene. The effect 
of the the electric field on the electronic system can be 
understood in terms of transfer of spectral weight, as we 
shall show in the next section. As the electric field is 
switched on, the spectral weight is reorganized, although 
in a way that still preserves the conservation rule im- 
posed by the /-sum rule, that was derived and analyzed 
in the context of low-energy graphene by Sabio et a/.—. 

As a last remark, since we are dealing with a system 
in which a polarized ac electric field is already present, 
one has to wonder about the possible influence of the 
polarization of the external probe to which the system 
responds. In fact, since we are analyzing the dynamical 
polarizability, which arises from the coupling between the 
electronic density and the potential induced by the exter- 
nal probe, the response of the system in the linear regime 
is insensitive to the polarization of the probe field. In or- 
der to see a response that depends on the polarization 
of the probe, we would have to analyze the response of 
the electronic current, which couples to the electric field, 
and whose linear response function is the conductivity. 
Notice, however, that the response function itself will 
not be altered due to this polarization, since in linear 
response it only depends on the properties of the sys- 
tem in the absence of the external probe by virtue of the 
fluctuation-dissipation theorem. 



Analytic approximations for the dynamical 
polarizability 



We now evaluate the dynamical polarizability (1351) us- 
ing the analytic approximations developed in Sect. Ill Bl 
We first consider the imaginary part of (l3"5"j) . This yields 
the response of the non-interacting system to an exter- 
nal probe of energy ui and momentum q, and is a build- 
ing block for the Random Phase Approximation (RPA). 
RPA is known to work well for doped graphene j 3 ^ where 
Landau's theory of the Fermi liquid provides a good de- 
scription of the low energy excitations. In the case of 
undoped graphene, the issue is more complicated due to 
the lack of screening near the Dirac points In our case, 
due to the opening of a gap at k = 0, we expect RPA 
to be sufficient to describe the main features of the re- 
sponse of the system once electron-electron interactions 
are taken into account. 



1. Fo -approximation 

In the i^approxirnation a gap opens up at zero mo- 
mentum, so it can be used for a qualitative description 
of the response to an external ac probe. We set T = 0, 
and for the moment we will restrict ourselves to undoped 
graphene. In this case, the dynamical polarization in- 
volves only the term that accounts for transitions be- 
tween the Floquet bands (0,—) and (0,+), i.e., transi- 
tions in which the number of photons is conserved: 



II Fo (q,u;) = 2^ + , - 

■ 1.1 C 1 -4- £ 



- \|2 



(37) 



c k+q : ' c k,0 



where |Xko) are the eigenvectors from the Fq- 
approximation, see Eq. We find 



Imn Fo (q,w) 
x (1 + 



1 



2„2 



4 ^ - F*q* 
A 2 



i 1 + ^-F v ) 6 " FW " ' (38) 



with A = V 'AA 2 + f2 2 — il being the gap opened at the 
first resonance. This is actually the dynamical polariz- 
ability for a "gapped graphene"; now the gap is due to 
the presence of the circularly polarized ac field. Gapped 
graphene has been studied extensively by Pyatkovskiy, 39 
who also derived analytical expressions for the real part of 
the polarization, for both doped and undoped graphene. 
We show the imaginary part of the polarizability in 
the ^"approximation in Fig. |4] For a given momen- 
tum of the external probe, the energy threshold re- 
quired to produce single particle excitations is increased 
due to the existen ce of the gap, being located now at 
oj = \J F^q 2 + A 2 . This yields a rearrangement of the 
spectral weight of the excitations, which might allow for 
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the existence of more complex excitations in the spec- 
trum of the interacting system. We investigate this ques- 
tion within RPA. The polarizability in the RPA is 



n 



RPA 



n (g,cj) 

1 - v q U {q., u) ' 



(39) 



where the denominator is the dielectric function in 
RPA with v q being the 2D unscreened Coulomb po- 
tential. In order to have long-lived collective excita- 
tions, i.e. plasmons, the dielectric function must van- 
ish at certain points u) p (q), which leads to the con- 
ditions WqRclIo(q, uS) = 1 and ImiIo(q, uS) = 0. In 
Fig. [5] the imaginary and real part of the polarizabil- 
ity in the RPA are plotted, where we use in Eq. (|3"9")l 
n (q, uj) = np (q, w) (Eq. (|37|)). It can be seen that 
the divergence in the threshold of excitations found for 
the non-interacting polarizability has disappeared, a fea- 
ture that is also observed in graphene in the absence 
of ac fields. However, the real part of the polarizabil- 
ity does not develop a resonance, leading to the ab- 
sence of plasmons, at least within the RPA of the Fq- 
approximation. In order to test if this is still true when 
higher order photon resonances are included, we analyze 
the Fi-approximation in the next section. 
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FIG. 4: Imaginary part of the polarizability in the Fo- 
approximation as a function of uj and for different values of q, 
Eq. ((38)) . For uj < ^/F ( 2 q 2 + A 2 , no excitations are possible, 
as compared to free graphene, where no excitations are pos- 
sible for uj < q. Note that only for small q = 0.05, the effect 
of the gap is visible as a shift of the divergence away from 
ui — q. Parameters: A = 0.1, Q = 1. 



Before that, let us analyze the case of doped graphene. 
As already shown in Ref. [39( for the case of gapped 
graphene, the plasmon already present in the system 
without ac fields is still robust once a weak field is in- 
troduced. The main effect of the external ac field is to 
modify the plasmon dispersion: 
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FIG. 5: Imaginary and real part of the polarizability in the 
Fo-approximation as a function of uj for different q in the 
RPA. The divergence in the imaginary part of the polariz- 
ability has disappeared, but the real part has not developed 
a resonance, which would be the signature of the existence of 
collective electronic excitations. Parameters: A = 0.1, = 1. 



^, = ,/«^(l 

(q) « y/q 

plasmon frequency is diminished due to the effect of the 
external ac field, since F < 1 and 1 — A 2 /(4^i 2 ) < 1. The 
correction coming from the factor Fq is essentially due to 
the renormalization of the Fermi velocity due to the ac 
field. The second correction depends on the relation of 
the chemical potential to the gap at zero momentum, and 
it is maximal for a chemical potential below A/ 2, where 
the plasmon is completely suppressed since no electrons 
are populating the upper Dirac cone. For chemical poten- 
tials above this value, the correction tends to be smaller, 
being almost negligible for /i ^ A/2. We point out that 
similar results hold in a quite different context, that of 
graphene anti-dot lattices ) 40 ' 41 where it was found that 
in the limit of low doping, gapped graphene models re- 
produce very well the plasmon dispersion of the anti-dot 
lattice. 

In short, the results from the ^"approximation yield 
a similar picture to that of graphene in the absence of 
an external field, and the main effect of the ac field 
is a renormalization of the single and many-particle 
spectrum, with a shift of the threshold for excitations 
due to the gap at zero momentum. 



A^-), but not the dependence on 
F t„\ „ rr. The 



2. F 1 -approximation 



'gNsNyFoq/i 



A 2 



l -T~S > ( 4 °) 



where g — e 2 /vpeoh is the fine structure constant of 
graphene. The correction affects the plasmon frequency 



The F\ -approximation accounts for non-zero-photon 
processes, neglected in the Fo-approximation. It is not 
fully analytically tractable in the calculation of many ob- 
servables. In the F\ -approximation using Eqs. (|25l) - (|2"5)l 
the polarizability for undoped graphene at T = be- 
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comes 



n Fl (q,w) = 



K£k+q,ol£k : o)| 2 

w - e k+ q ,o + e k~,o + iV 



+K^ +q ,-ilCo> 



+K^ +q ,-2l€k,o> 



1 



1 



^ - e k+ q ,-i + e k,o +*V e k+q 4 + e k + n 1/ 

1 1 



^ - £ k+q,-2 + £ k,0 + l7 l £ k+q,2 + £ k,0 + *V, 



(41) 



Here, there are three different contributions to the polar- 
izability, coming from the Floquet bands I = 0, I = ±1 
and I = ±2, i.e., from excitations that involve the ex- 
change of up to two photons from the external field. How- 
ever, for the electric fields in which this approximation 
holds, the contribution from the I — ±2 components is 
essentially negligible, and therefore only zero and one- 
photon processes will be considered. In what follows we 
evaluate ([41 [I numerically, first integrating the imaginary 
part and then computing the real part via the Kramers 
Kronig relations. 

Figure [6] shows the imaginary part of the polarizabil- 
ity Hf 1 for fixed q as a function of w. In the upper 
panel, the components I = and I = ±1 and their sum 
are shown for q = 0.1, in order to illustrate where its 
structure comes from. The two lower panels represent 
the total polarizability for two different wavevectors q, 
divided into two regions for better visibility of the dif- 
ferent features. Several new features emerge from the 
Fi -approximation. As shown in Fig. [BJ at the level of 
zero-photon processes, there is a gap at zero momentum, 
already captured in the Fq— approximation. In addition, 
gaps appear at higher momenta, where the first anticross- 
ing of Floquet sidebands occurs (see Fig. For a suffi- 
ciently small q this second gap translates into two small 
gaps in the single particle excitation spectrum around 
w « 1, which are eventually closed for higher momenta, 
as shown in Fig. [5] (lower panel, right plot). The first of 
those gaps, for 0.9 < w < 1 , is due to the fact that for 
electrons from the lower cone of graphene no states are 
available in the upper band for those values of q and ui 
due to the anticrossing of Floquet sidebands. The second 
one, at 1 < w < 1.1, is due to the lack of states in the 
lower cone in the region where this anticrossing occurs 
with lower Floquet sidebands. 

The most important new features of the response of 
the system, however, come from the contribution of one- 
photon processes, in which transitions from the I = to 
the I = ±1 sidebands are taken into acc ount. New single 
particle excitations appear below u = yjF^q 2 + A 2 , leav- 
ing only a small region of energies where no excitations 
are found, a region which again is closed for sufficiently 
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FIG. 6: Imaginary part of the polarizability IIfi as a function 
of cj. Upper panel: q = 0.1; the components I — 0, ±1 of the 
polarizability and their sum are shown, see Eq. [?T] Lower 
panels: Imaginary part of the total polarizability as a function 
of ui for two different q: q\ = 0.05, q-2 = 0.1. Left plot for 
cj < 0.4, right plot for u) > 0.7. Parameters: A = O.f , Q = 1. 



large q (dashed line for q = 0.1). 

One-photon processes introduce new single particle ex- 
citations into the response of the system, and we next 
examine the effect of these processes on the collective ex- 
citations of the system. The RPA polarizability Hp 1 ; rpa 
is shown in Fig. [7] for different values of the external mo- 
mentum q. One-photon processes have an important 
effect on the response of the interacting system, allow- 
ing for the existence of collective excitations for small 
enough momentum, see curves for qi = 0.035, 92 = 0.05 
in Fig. [7] For those, the plasmon conditions arc fulfilled, 
which is reflected in the development of a resonance in 
the real part of the RPA polarizability. For an external 
momentum of q — 0.05 e.g., the resonance is located at 
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FIG. 7: Imaginary and real part of the polarizability IIpj 
in the RPA approximation. Notice the resonance in the real 
part for small momenta q\ = 0.035 and qi = 0.05 (lower 
panel), where no excitations in the imaginary part exist (up- 
per panel), pointing to the existence of collective excitations. 
Parameters: A — 0.1, Q = 1. 



lj ~ 0.021, which is already in the region where the imag- 
inary part of the polarizability is zero, allowing for an 
undamped plasmon. It is important to remark that this 
plasmon becomes unstable in two different scenarios, (i) 
For large enough momentum of the external probe, where 
the resonance is weakened and it occurs in a region where 
single particle excitations exist, so the plasmon can de- 
cay into those excitations, see q 3 — 0.15 in Fig. [7] (ii) 
When two-photon processes are considered, there is no 
longer a region where the imaginary part of the polariz- 
ability is zero. For weak fields, however, these processes 
are negligible and their effect on the plasmon should also 
essentially be irrelevant. However, this suggests that as 
we increase the intensity of the electric field, and higher 
order photon processes are important, there is no region 
of momenta q in which the plasmon is stable. 

For doped graphene, the -Fi -approximation introduces 
similar features as those described for undoped graphene, 
see Fig. [5J The effect of the anticrossing of Floquet side- 
bands is to induce gaps in the response of the system 
for lj 1, and processes including the exchange of one 
photon give rise to new excitations for small energies. In 
order to quantify the effect of these processes, in this fig- 
ure the polarizability is compared to the one for doped 
graphene, using the Fq— approximation, where only zero- 
photon processes are considered. The RPA response of 
the interacting system is shown for a couple of represen- 
tative external momenta in Fig. [9j As it happened in the 
undoped case, for small external momentum (q = 0.05 in 
Fig. in]) there is a resonance in the real part, signaling the 
existence of a plasmon, which again has a renormalized 
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FIG. 8: Imaginary part of the polarizability for doped 
graphene as a function of lj for q — 0.05. Both the Fq- (solid 
red) and the f\-approximation (dashed dark red) are shown. 
In the lower panels, the upper plot is split in two parts in or- 
der to better visualize the different regions of lj. Parameters: 
A = 0.1, fi = 1, fi = 0.2. 



dispersion relation due to the effect of the external ac 
field. However, non-zero-photon processes are respon- 
sible again for the appearance of low-energy excitations 
that tend to make the plasmon unstable for large enough 
momenta q {q = 0.15 in Fig. [9]) and for larger intensities 
of the field, as discussed for the undoped case. These 
momenta q for which plasmons become unstable are still 
lower than those for which the plasmon is damped in 
graphene without ac field. 

Summarizing, the inclusion of non-zero-photon pro- 
cesses is crucial in order to capture the physics of the 
response of graphene to an external probe in the pres- 
ence of a weak ac field. This is due to the appearance 
of excitations in the low energy spectrum of the system, 
not included in the ^"Approximation, that allow for the 
existence of collective excitations in undoped graphene, 
but make those plasmons unstable for smaller momenta 
than their counterparts in graphene with no ac fields. 



IV. CONCLUSIONS 

In this work we have analyzed the properties of 
graphene under an external circularly polarized ac field 
in the weak field regime. We have developed analyt- 
ical approximations to the Hamiltonian, the so called 
Fq (see section III B 1|) and F\ -approximations (see sec- 
tion QTB2]), that allow a certain analytical tractability of 
many relevant objects. The i<o _a PProximation includes 
only zero-photon excitations in the system, and is useful 
to calculate certain observables in the low energy sec- 
tor. The F\ -approximation includes higher order photon 
processes, allowing for the analysis of a wider range of 
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FIG. 9: Polarizability in the RPA approximation for doped 
graphene in the Fi -approximation. Upper panel: imaginary 
part. Lower panel: real part. The results are shown for two 
different momenta q — 0.05 and q = 0.15. In both figures the 
results are compared with the Fo-approximation. Notice the 
existence, in both approximations, of the resonance in the 
real part of the polarizability, signalizing the existence of a 
collective excitation. For large momentum q = 0.15, however, 
the imaginary part shows no gap and therefore the plasmon 
can decay into single particle excitations. Parameters: A — 
0.1, il = l,n = 0.2. 



observables and a larger energy sector. However, it re- 
quires in many cases numerical calculations to extract 
the observables. 

Special emphasis has been put on the calculation of 
the polarizability of the system, which can be used to 
analyze the spectra of single and many-particle excita- 
tions of the system. We have derived a general expres- 
sion for the polarizability of graphene in the presence of 
an ac electric field, which we have analyzed in the con- 



text of the Fq and F\ -approximations. While the for- 
mer allows for analytical expressions, and captures well 
the effect of the zero-momentum gap in the system, it 
misses the non-zero-photon processes, that are captured 
by the Fi-approximation, and in turn are responsible for 
the emergence of collective excitations even for undoped 
graphene, as far as the Random Phase Approximation 
remains valid. However, it also points out that these 
collective excitations are less stable when compared to 
graphene with no external ac field: for large enough ex- 
ternal momenta and ac field intensities, these excitations 
become damped and acquire a finite lifetime. 

We have shown that circularly polarized ac fields can 
be used to modify the properties of graphene in several 
ways: (i) They open up gaps at zero momentum that can 
be exploited in practical applications, (ii) They permit 
the existence of plasmons (in both undoped and doped 
graphene), (ii) The plasmon frequency is tunable with the 
external field, and, finally, (iv) For large enough fields the 
plasmons become unstable. Moreover, we have developed 
and tested analytical tools to analyze theoretically the 
behavior of graphene in the presence of ac electric fields, 
which should be useful in future works in this field. 
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Appendix A: Derivation of the polarizability for circularly polarized field 

The derivation of the formula for the dynamical polarizability follows the lines of its counterpart in the 2DEG. 28 
The wavefunction for graphene under a periodic driving can be written by use of the Floquet theorem as 

T/Wr, t) = -ie^e-^'^W, (Al) 
V2 

where is the quasienergy and 0k,<r(*) are the Floquet states which fulfill the time-periodicity of the driving field, 
and we have chosen the solution corresponding to the First Brillouin zone. After applying a weak probe potential, 
these wavefunctions are not any more eigenfunctions of the full Hamiltonian, but we can use them as a basis to write 
the new wavefunction: 

*k,„(r,t) = ^a k >/(t)Vk', CT '(r,i) ( A2 ) 

kV 
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Inserting this into the Schrodingcr equation for the Hamiltonian Ho(t) + Hi(t), where Ho(t) is the Hamiltonian for the 
periodically driven graphene, and H\{t) = V(r,t) represents the weak probe potential, we are left with a differential 
equation for the coefficients a^ a {t): 



i Y " k '> CT ' (*)V>k', CT ' ( r ' *) = Y ak '' a ' W^( r > ^V'k'.a' (r, t) 
kV kV 

We can now project this equation into a state Vk",o-"i yielding 

a k »,„»(*) - -i^ak-, ff '(i)e i(6k "-"-^.'' )1 ft v ,(i)^,,r'(i)7(k" -k',t), 

kV 



(A3) 



(A4) 



where V(k" — k',t) is the projection of the probe potential into the states k' and k". Now we can expand this 
equation in a power series of the external potential, and keeping only the first order we are left with 



This equation can be simplified by Fourier transforming it, yielding 
o k //, ff »(t) - y 7^(k" - k, j)e- iut e i ^"-"- ek ^ t e r > t Y 



±n f .a* in. a* , in'.b* ±n.b* 



oj - (ri - n)Q - (ek",cr" - ek,<r) + *»7 



(A5) 



(A6) 



In order to get the response of the system to the external probe in linear response, we write down the expression 
of the induced charge density: 



pj£(r, t) = *£ iCT (r, 0*k, CT (r, t) - ^(v, t)^ a (r, t) 

= Y a k V (*)^k',<7' ( r . ^V'k.a (r, t) + o k ', CT ' (t)^ kjCT ( r , i)V"k', CT ' (r, t) 



kV 



and insert the result obtained for a^ t<T {t). After some algebra we arrive at 

r, J o*> k 



(A7) 



(A8) 



where we have introduced the short notation 



nn' mm' 



1 

2 

if 
2 



„i(n' -n)nt Aim' -m)nt ( jn' ,a* <n,a . in',b* <n,b\ ( <m' ,a* ,m,a 
e e ^k+q.CT'^k.ff "+" V > k+q,(T' < Pk,(T y / ^k,<7 ^k+q.cr' 

w - (n' - ra)fi - (ek+ q ,<T' - ek )C r) + ^ 

i(n'-n)nt e -i{m'-m)nt ( A n ' ^ a (h n ' a * 
\ ' k — q.cr' ' k,c 



+ 



im f ,b* ±rn,b 



n',6 /n.6*\ / ±m'a im.a* . / m 
k-q^'^a ) (<Pk,<7 •Pk-q.ff' + <?k 



cr ^k— q,<ry 



-w - (n' - n)fi - (ek-q,<r' - £k, CT ) - i»? 
By comparing with the Poisson equation 

p ind (r,t)=Y [ ^y ind (q, w )e-^ t e lqr |- 
q 

we see that the induced potential must fulfill 

T/ ind (q,c) = ljl/^(q,c)^^/ k , ff J-k,^. 
q era' k 

We sum on both sides y ext (q, w) and get 



• (A9) 



(A10) 



(All) 



(A12) 
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where the dielectric function is given by 



e(q,") = 7-4^ — — = • ( A13 ) 

In the RPA approximation we obtain therefore 

e(q, w) RPA = 1 j X! X! fKe^k,*,*' ■ (A14) 

erer' k 

After substituting the expression for J-k.o-.u' and some straightforward manipulations, we arrive at our desired result 
for the dynamical polarizability: 



, , w - ek+q, CT ',; + Ck.tr.o + l V 

oa' k / ' ' 



Notice that now we have simplified the expression by writing it as the scalar product between different Floquet 
sidebands by using Eqs. ([5]) and ([7|). 
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